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Abstract 

In this paper we use quantum magnetohydrodynamic (MHD) as well as magneto- 

£3 ■ 

hydrostatic (MHS) models for a zero-temperature Fermi-Dirac plasma to show the 
fundamental role of Landau orbital ferromagnetism (LOFER) on the magnetohydro- 
static stability of compact stars. It is revealed that the generalized flux-conserved 
equation of state of form B = j3p 2s ^ only with conditions < s < 1 and < (3 < \p2m can 
leads to a stable compact stellar configuration. The distinct critical value f3 cr = \^2tt is 
shown to affect the magnetohydrostatic stability of the LOFER (s = 1) state and the 
magnetic field strength limit on the compact stellar configuration. Furthermore, the 
value of the parameter f3 is remarked to fundamentally alter the Chandrasekhar mass- 
radius relation and the known mass-limit on white dwarfs when the star is in LOFER 
state. Current findings can help to understand the role of flux-frozen ferromagnetism 
and its fundamental role on hydrostatic stability of relativistically degenerate super- 
dense plasmas such as white dwarfs. 

PACS numbers: 52.30.Ex, 52.35.-g, 52.35.Fp, 52.35.Mw 
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I. INTRODUCTION 

Since the pioneering discovery by Chandrasekhar in 1939 [l| concerning the 
mass-limit on compact stars and the hydrostatic stability mechanisms in such 
stars due to the relativistic degeneracy of electrons, there has been a growing 



* degenerate ionized 
6j. It is well known 



interest towards the study of hydrodynamic properties o 
matter, the so-called zero-temperature quantum plasmas {2- 
that the matter under compression exerts enormous pressure called the degen- 
eracy pressure due to the Pauli exclusion principle when the interparticle dis- 
tances are lowered to become comparable to the de Broglie thermal-wavelength 
X D = h/(2Trm e k B T) 1 ' 2 [7). Quant um peculiarities such as quantum tunneling, 
quantized Hall effect, magnetic quantization etc. ubiquitously appear as the 
degeneracy limit is reached. Such peculiar features prove to be of fundamental 
significance and applications in ordinary metallic and semiconductor materi- 
als. Many recent investigations based on the quantum hydrodynamics (QHD) 



and quantum magnetohydrodynamics (QMHD) models |8l-4l7l] indicate that the 
incorporation of the quantum electron-tunneling and degeneracy pressure can 
lead to quite different nonlinear dynamic effects in plasmas. It has also been 
remarked that the relativistic degeneracy caused by large-scale gravitational 
forces in stars which causes the gravitational collapse in stellar objects jisl . [nJ ] 



may also lead to distinctive nonlinear hydrodynamic features [2CH22J] due to the 
change in the thermodynamical quantities in the Fermi-Dirac statistics 

Among the greatest challenges today is the problems associated with the ori- 
gin of strong magnetic fields present in many compact astrophysical entities such 
as white dwarfs, pulsars, neutron stars, etc. and its formidable role on the stellar 
chain of evolution. It is also believed that the magnetic field has a fundamental 
role in the formation and the dynamical processes in astrophysical environments 
(24I ]. There has been extensive past studies on the thermodynamical behavior 



of degenerated electron gas under arbitrarily high magnetic field [25-32]. Such 
investigations have revealed that the high magnetic field can lead to the anoma- 
lous quantization and spiky features in the electronic density of states (DoS) 
affecting all the thermodynamical properties of the Fermi-Dirac gas. It was sug- 
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gested that under such a quantizing field the magnetic transverse collapse of the 



gas is possible 
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where the Fermi-Dirac gas may become one-dimensional. 



A review of the current findings on the properties of matter in strong magnetic 



361 ]. Recent studies based 



field has been reviewed in some recent literature [35, 
on the QMHD including the magnetization effect confirm the significant role of 
the electron spin-orbit magnetization effects on the nonlinear properties of de- 
generated quantum plasmas [37H47I]. Particularly, a more recent study remarks 
distinctive paramagnetic nonlinear features of a Fermi-Dirac plasma due to the 



relativistic electron degeneracy 



48|]. 



Another outstanding feature of a degenerated plasma under a strong magnetic 
field is that a ferromagnetic solution called the Landau orbital ferromagnetism 
(LOFER) is possible [491-1541] which may account for the large magnetic fields 
(as high as 10 8 G) estimated for some astrophysical compact objects. In the 
present investigation we use both magnetohydrodynamic and magnetohydro- 
static (MHS) models to explore the role of Landau orbital ferromagnetism on 
compact stellar characteristics such as hydrodynamic quantum collapse, hydro- 
static stability and the mass-radius relation comparing the results with that of 
previous ones. The presentation of the paper is as follows. The QMHD model 
including Bohm potential and the spin-orbit magnetization effects is introduced 
in Sec. HI] and the possible regimes for a quantum collapse is explored in Sec. 
IIIIl The hydrostatic stability of a LOFER state is investigated and the Chan- 
drasekhar mass-limit is calculated based on the generalized LOFER equation of 
state in Sec. IIVI Finally, a summary is given in Sec. |V] 



II. MAGNETOHYDRODYNAMIC MODEL AND DEGENERACY COLLAPSE 

In this section we show that for a flux-conserved degenerate plasma the total pressure can 
vanish in the transverse direction making the MHD wave unstable due to the effect of plasma 
magnetization. In this case a transverse collapse may occur without the need for presence 
of the gravitational force. Let us consider the magnetohydrodynamic (MHD) equations for 
a completely degenerate dense singly-ionized and quasineutral (rij ~ n e = n) helium 
plasma with the center of mass density p = m^ni + m e n e ~ m^n = 2m p n (m p is the 
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proton mass). We have for the continuity equation 
dp 



Of 



+ V • (pu) = 0, 



(1) 



where, u = (m^Ui + m e riiU e ) / p is the center of mass speed of plasma. Therefore, the 
generalized momentum equation including the quantum degeneracy pressure, dipole force, 
magnetization and electron nonlocality effects can be written in cgs units in the following 
form 



p| = (B-V)(B-M)-v(y-M-B 



VP d +^- V- v 



(2) 



2m e m; y/p 

in which is the electron degeneracy pressure and we have neglected the Bohm 
force on ions and the pressure due to them. The magnetization, M and the in- 
duced field, B, are related through; H = B — 47rM(B), with H being the magnetic 
field due to physical currents. The metastable Landau orbital ferromagnetism 
(LOFER) for the flux-conserved degenerate plasma model is given by H = 
which leads to B = 47rM(B) [55[. On the other hand, it has been shown that 
the LOFER condition for a magnetized degenerate electron-gas leads to the 
field/density equation of state (EoS) of form; B{r) = a(pe/ p e ) 2//3 |50| . where, 
P6 = p/10 6 and jj, e is the number of nucleon per electron (in this calculation we 
use p e = 2 for helium). The parameter a is the normalizing factor for magnetic 
field and will be found to be of the order 10 s for the white-dwarf mass-density 
ranges (p 6 ~ 1). However, this parameter is known to relate to some other pa- 
rameters such as the plasma temperature, electron exchange interactions etc. 
1 531]. There has been many reports of compact star with strong internal or ex- 
ternal magnetic fields |56l-459l]. On the other hand, many reports confirm the role 
of flux-conservation on star formation | 63_] with the similar density dependence 
of magnetic field. In 1964 a theory based on flux conservation was suggested 
simultaneously by Ginzburg [6o| and Woltjer [61] to explain the presence of 
intense magnetic fields in some young compact stars born in supernova explo- 
sions | 62| . As it will be apparent later in discussion the parameters s and a 
in EoS of form B = ct(p 6 /p e ) 2s / 3 are central to the stability criteria of compact 
stellar configurations. There are some theoretical discrepancies on the values 
of these parameters. For instance, some calculations show that [54] the value 



of s in LOFER EoS, should be 4/3 rather than 2/3 calculated for this parame- 



ter previously [501 ]. Also, the calculations of the parameter a leads to different 
values due to the oscillatory nature of the spin-orbit magnetization elements. 
However, in this calculation we introduce a more general field/density EoS as 
B = a(ps/ p e ) 2s l 3 to show that only the restricted values of < s < 1 is consistent 
with the magnetohydrostatic (MHS) stability of compact stars. Also, we will 
find an upper limit on the value of a (or the strength of the magnetic field) for 
the known LOFER state (s = 1) of stellar configuration. 

In a strongly magnetized Fermi-Dirac electron gas the equation of state (EoS) is quantized 



64] 



and we may write for the electron number- density in terms of Hurwitz zeta functions 
n e {x )1 )=n c {2 1 fl 2 H_ 1/2 (g) , 

H g (q) = h(z, {q}) - h(z, q + 1)- \q~\ 

oo 

H z ><l) = E {n + qy z - 

n=0 

where h(z, {q}) is the Hurwitz zeta-function of order z with the fractional part of q as 
argument and P e » denote the degeneracy pressure parallel to the magnetic field. Note also 
that, n c = mlc 3 /2ii 2 h 3 , 7 = B /B c with B c = m 2 e (? /eh ~ 4.41 x 10 13 G being the fractional 
critical-field parameter, ep e = y/1 + x 2 = EF e /m e c 2 is the normalized Fermi-energy and 
x = pFe/meC is the normalized Fermi-momentum the so-called relativity parameter. It is 
evident that, the electron degeneracy pressure is field dependent on both magnitude and 
direction so that P e ± ^ P e ». It has been shown that in such plasmas two distinct quantum 
and classical degeneracy regimes based on the parameter x 2 /2 7 can be defined [65]. In our 
case the magnetic fields of interest are of the order 10 8 G leading to the classical regime 
x 2 /2'j ^> 1 for the white dwarf star density ranges. Hence, one may assuredly use the 
Chandrasekhar classical EoS for the degeneracy pressure as [l| 
4 5 

p d {x) = jx (2x 2 - 3) VTT^ 2 + 3 In [x + Vl + x^\ } , (4) 

where, the well-known relativity parameter is x = j9F e /m e c = (h/m e c)(3n/8Ti) 1 ^ 3 = 



(n/rio) 1 / 3 (no = 8irmlc 3 /3h 3 ~ 5.9 x 10 29 /cm 3 ) |23l ]. with p Fe being the relativistic 
Fermi momentum. We may also write the relativity parameter in terms of 
the plasma mass-density p = 2m p n and p = 2m p n = 2m p / 2>it 2 X 3 c ~ 2 x 10 6 gr/cm 3 



as x = (n/n o y/ 3 = (p/po) 1 ^ 3 with A c = h/m e c ~ 3.863 x W~ n cm being the scaled 
electron Compton-wavelength. Note that in the forthcoming algebra we will use 
the normalized mass-density p = p/po dropping the bar notation for simplicity, 
hence, x = p 1 / 3 . Therefore, the normalized momentum equation, Eq. ([2]) may be cast in 
terms of the effective potentials [66J as 



p— = -V 



ph 2 „A*/p 



(5) 



2m e mj */p 

where \l/ m and ^ are the corresponding flux-conserved ma gnet ic and electron degener- 
acy effective potentials (the normalization is given elsewhere 66]). Therefore, the inward 
magnetic-force, F m = V^ m = — sCl(5 2 Vx^ s ~^ /2tt(As — 3) (we have made use of a = /3C S 
with C s = c^Jm e /2m p being the quantum sound speed) and the outward degeneracy force, 
Frf = V 1 !^ = C 2 Vyl + x 2 oppose each other. It has been shown that on some conditions 
the total force is canceled giving rise to the quantum collapse for LOFER case 66j. However, 
for the general EoS considered here there exists a surface for which a quantum collapse can 
occur which is shown in Fig. 1 for the parameter range of 3/4 < s < 1 . It is clearly observed 
that for a quantum degeneracy collapse to occur in the LOFER state (s = 1) the value of 
should be above a critical value (3 cr = \/2tt 66]. It is also remarked that quantum degener- 
acy collapse is not possible for s < 3/4. In the proceeding section we use the MHS stability 
of a magnetized degenerate plasma model to find some limits on the value of parameters s 
and f3 which is also a limit on the magnetic field of flux-conserved compact star. 



III. HYDROSTATIC STABILITY IN FLUX-FROZEN STELLAR PLASMAS 

Using a classical treatment, in this section, we consider the MHS equilibrium for ho- 
mogenous spherical magnetized plasma. We also ignore the quantum tunneling effect on 
electrons due to large size of system compared to the interparticle distances. To begin with, 
let us consider a spherical shell-element of mass-density p which is normalized to po defined 
in previous section. The equilibrium for the shell-element is defined as Fi n = F out , where, 
the inward force consist of the gravity and the magnetic force and the outwards force is the 
electron degeneracy force. In a simpler argument we may write 

dP tot Gp(r)M(r) 



dr r 



■ Ptot — Pd + Pm-, (6) 
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where, G is the gravitational constant and M(r) is the mass of plasma within the radius r, 
defined as 
dM{r) 



dr 



Attt p(r) 



(7) 



where, again, p is the normalized local density. Thus, the above definitions lead to the MHS 
stability condition of the form 

I d ( r 2 dPt t(p) dp(x) dx{r 



+ p(r) 



0. 



(8) 



4nr 2 G dr \p(r) dp dx dr 
with the parameter x = p 1 ^ 3 being the Chandrasekhar relativity parameter. From the 
definitions for the degeneracy and magnetic pressures one obtains 



dPdjp) 
dp 



C 2 S 



x 



6 VT+ 



xs 



dPrn(p) 

dp 



12tt 



-x 2s ~\ (9) 



where, we have used a general form of B 
state previously defined. Hence, we have 



ap 



2s /3 



[s is real) for the magnetic equation of 



dPtot(p) C 2 



dp 



6 



x 



_VT+ 



2tt 



-X 



2s-l 



Thus, in terms of the relativity parameter, we rewrite Eq. (|8]) as 



3iiX 2 n 2 h d 
16r 2 dr 



x{r) 



l + x(r) 



2tt 



-x(r 



,2(a-l) 



dx{r) 
dr 



+ x{ry 



(10) 



(ii) 



where, ny L = ^hc/G/2m p ~ 1.3 x 10 19 is the dimensionless hierarchy-number defined based 
on three fundamental constants, namely, the scaled Plank- const ant, h, the speed of light in 
vacuum, c and the gravitational constant, G. It is clearly evident that, for the value of s = 1 
the magnetic field pressure gradient term in Eq. ( TTTT) becomes independent of the relativity 
parameter. 

Therefore, each stable model is obtained by integrating Eq. (11 ip from the center (r = 0) 
outwards until x(r) vanishes for some r = r s which would be the surface (of the star). This 
is done by employing initial conditions x(r = 0) = x c and dx/dr(r = 0) = so that for every 
given value of the relativity parameter at the center (x c ) we find a distinct solution. However, 
it is convenient to put the Eq. (TTTT) in a more friendly shape by introducing dimensionless 
parameters y = x/x c and r\ — r/r c with r c = v / 37rn/jA c /4x c so that at the center (77 = 0) we 
will have y — 1. Therefore, the Eq. (1TT|) in dimensionless form reads as 



1 d 

rf drj 



\/ 1 + x ly(vY 



s/3 2 
2tt 



-x 



2{s-l) 



y(v) 



2( S -1) 



dy(y) 
dr] 



+ y(r ] ) 3 = 0, 



(12) 



It is noticed that at the very high central density (x c — > oo) Eq. ( IT2|) reduces to the famous 
Lane-Emden equation of index 3, which for /3 = (unmagnetized case) leads to the well- 
known Chandrasekhar mass-limit of Men — 1.43Msf (Ms being the mass of the sun). A 
standard integration algorithm such as the Runge-Kutta fourth-order (as I have carried out 
the calculations with mathematica software) may be used to find, for instance, the value of 
7] = 7] s at which y vanishes (the surface of star) from which the radius and the mass of the 
stable configuration (star) can be calculated via the following relations 



3tt 



>ls 



r c r} s , M = An J p(r)r dr = - m p n h J y(rj) rj drj. (13) 



o 



It is interesting, however, to note the unique scaling, i.e. rihX c ~ 0.787Re (Re being the 
earth's radius) and 2m p n\ ~ 1.8A9Ms (Ms being the sun's mass). These relations for (3 = 
(unmagnetized) case, following the Chandrasekhar's pioneering work, has been reviewed 
in Ref. 67|. The main goal of next section is to explore the effect of Landau orbital 
ferromagnetism on magnetohydrostatic stability and consequences on the Chandrasekhar 
mass-limit in completely degenerate magnetized plasmas. 



IV. THE CHANDRASEKHAR LIMIT FOR FLUX-FROZEN COMPACT STAR 

Equation (TlT]) can be solved for given values of 0, s and x c . The solution y(rj) for 
the simplest case of /3 = (or s = 0) which corresponds to the unmagnetized (uniformly 
magnetized) plasma case is shown in Fig. 2 for different values of central fractional density 
x c . The values of t]q where y(t]o) = indicate the surface of plasma from which the radius 
and mass of the configuration can be calculated. It is observed that as the value of x c 
is increased r] also increases and in the limit it approaches a liming value r]ch called the 
Chandrasekhar value. Figure 3 shows the famous Chandrasekhar mass-radius plot (in terms 
of sun's mass and earth's radius) and the corresponding Chandrasekhar mass-limit Men — 
I.43M5 for unmagnetized case. Following is the evaluation of the problem for different values 
of parameters s and 0. 



A. Case s > 1 



For the case s > 1, it is observed (e.g. see Fig. 4 for s = 2) that, as one increases the 
value of x c to x c — 1, the value of rj increases and approaches to a liming value. However, 
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it is remarked that, there are no stable solutions for x c > 1 in this case. In other words, in 
this model, i.e. s > 1, no dense-centered stable configuration such as observed for a star is 
allowed. Therefore, this case can not correspond to a physical solution. 

B. Case s < 

On the other hand, for negative values of s, when x c is finite no 7/o-value can be obtained, 
hence, the plasma must be infinite. However, for x c = +oo the same Chandrasekhar-limit is 
obtained, regardless of the value of parameter /3, as it is evident from Eq. (jlip . For instance, 
Fig. 5 shows the solution y(rj) and the variation with respect to different values of x c for 
the case s = — 1, indicating that, the case s < is unphysical. 

C. Case < s < 1 

For the range of < s < 1 the plasma is stable for large values of x c . However, as 
x c — ^ it is observed from Eq. ( ITT]) that the integrand diverges and the plasma becomes 
unstable. Therefore, in this model there may be a radius-limit to the plasma evaluation of 
which is beyond the scope of this paper. Figure 6 for s = 1/2 shows the similar features as 
Fig. 2 as the central density of plasma is increased. Also, evaluation of Eq. ( TTTT) reveals 
that the same Chandrasekhar mass-limit is obtained and that in this model the mass-limit 
is invariant under the change in /3 parameter (e.g. see Fig. 7). Another important feature 
in this model is the existence of critical value for /3. Evaluation of the solutions to Eq. ( ffTj) 
for this case shows that stable configurations exist only below a critical value f3 cr = 2n. 

D. The LOFER Case s = 1 

In this model which corresponds to the equation of state of form B = (3p 2 ^ 3 and is known 
as LOFER state, unlike the previous case, the plasma is stable for the whole range of the 
parameter x c . However, the stability of the plasma is removed for values of (3 > (3 cr . It is 
also remarked from Fig. 8 that in this model the chandrasekhar-limit decreases as the value 
of is increased. From the previous discussion it may be concluded that the only stable and 
physical model corresponds to the cases s = 1 and (3 < 2n. The magnetogravity collapse 
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condition (0 < 2n) is to be compared with the one obtained for the quantum collapse 
mentioned above (0 > 2tt). Therefore, in a dense ferromagnetic plasma depending on the 
value of only one of possible collapses is possible of which the magnetogravity collapse can 
occur in hydrostatically stable configurations. It is noted that, although the LOFER-state 
is a metastable configuration and requires a sufficient condition to operate, however, once 
it occurs the quantum collapse will inevitable. The variation of Chandrasekhar limit with 
the LOFER parameter 0, shown in Fig. 9 confirms that the mass-limit decreases with the 
increase in the value of until it vanishes at critical value = \/2~k. Also, Fig. 10 compares 
the Chandrasekhr mass-radius curves for the tentative values of 2 /2n = 0,0.1,0.5 and 
s = 1. 

It is remarked that, the flux frozen assumption used in this investigation can lead to 
fundamental effects of the Chandrasekhar mass-radius relation for the case of s = 1 and the 
mass-limit on magnetic white dwarfs can be much lower than 1.43217. This stellar model is 
also consistent with supernovae type I theory of Hoyle and Fowler {(38 1 which requires stellar 
mass between 1.01 and 1.40 solar mass. It is obvious that for a complete treatment of stellar 
stability the effects such as Thomas-Fermi, Coulomb, electron exchange and ion correlation 



have to be considered 



69]. Also, the electron-capture (inverse beta-decay) has been shown 



to alter the known EoS and limits on the compact stars [70] investigations of which is out 
of the scope of the current study. 



V. CONCLUSION AND SUMMARY 



Using the quantum magnetohydrodynamics and magnetohydrostatic models we showed 
that in a relativistically degenerate magnetized stellar objects such as white dwarfs the 
magneto-gravitational collapse is possible. The collapse was shown to be consistent with 
Landau orbital ferromagnetism with magnetic equation of state B = 0p 2s ^ 3 on some con- 
ditions. It was also revealed that, a critical cr = y/2~Tr exists above which the plasma is 
hydrostatically unstable which comparing with the previous results indicates that in the 
LOFER plasma state quantum collapse and magnetogravity collapse can not take place 
at the same time. Furthermore, in a stable LOFER plasma it was shown that the Chan- 
drasekhar mass-limit may critically depend on the value of the and s-parameters. These 
findings reveal the inevitable role of magnetism on the evolution of dense stellar objects. 
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FIGURE CAPTIONS 



Figure- 1 

The surface in (3-x-s space showing where the quantum collapse is possible. 
Figure-2 

Variation of the value of r] with respect to change in the central density parameter x c for 
unmagnetized ((5 = 0) case. The thickness of the curves is used as a measure for the value 
of the varied parameter. 

Figure-3 

The famous Chandrasekhar mass-radius curve for the case of unmagnetized (ft = 0) 
plasma indicating the Chandrasekhar mass-limit of Mch — 1.43Mg. 

Figure-4 

Variation of the value of r] with respect to change in the central density parameter x c 
for magnetized case (B = /3p 2s ^) with (s = 2) showing magnetohydrostatic instability of 
this model for x c > 1. The thickness of the curves is used as a measure for the value of the 
varied parameter. 

Figure-5 

Variation of the value of rj with respect to change in the central density parameter x c 
for magnetized case (B = /3p 2s ^) with (s = —1) showing unphysical results for this model. 
The thickness of the curves is used as a measure for the value of the varied parameter. 

Figure-6 

Variation of the value of rjo with respect to change in the central density parameter x c 
for magnetized case (B = /3p 2s ^) with (s = 1/2) showing the magnetohydrostatic stability 
for the range of central plasma density parameter, x c ^ 0. The thickness of the curves is 
used as a measure for the value of the varied parameter. 

Figure-7 
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Variation of the value of r/o for fixed large central density parameter x c and varied mag- 
netic parameter, /3 2 < 2tt, for magnetized case (B = /3p 2s ^ 3 ) with (s = 1/2). The plot 
indicates that in this model the Chandrasekhar mass-limit does not depends on the value of 
the magnetic parameter, (3. The thickness of the curves is used as a measure for the value 
of the varied parameter. 

Figure-8 

Variation of the value of r] for fixed large central density parameter x c and varied mag- 
netic parameter, (3 2 < 2tt, for magnetized case (B = (3p 2s ^ 3 ) with (s = 1) showing the 
magnetohydrostatic stability for the whole range of central plasma density parameter, x c . 
The plot also shows that in this model the Chandrasekhar mass-limit depends strongly on 
the value of the magnetic parameter, j3. The thickness of the curves is used as a measure 
for the value of the varied parameter. 

Figure-9 

Variation of the value of the Chandrasekhar mass-limit with the value of the ferromagnetic 
parameter, j3 for the case of s = 1. The horizontal line indicates ordinary mass-limit 
M/M s ~ 1.43217. 

Figure- 10 

The Chandrasekhar mass-radius curves for ordinary LOFER (s = 1) with /3/V2tt = 
(rectangles), /3/y/2~7T = 0.1 (empty circles), and I3/^/2~k = 0.5 (filled circles) for the case 
s = 1. 
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